Density of states in an optical speckle potential 
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We study the single particle density of states of a one-dimensional speckle potential, which is 
correlated and non-Gaussian. We consider both the repulsive and the attractive cases. The system is 
controlled by a single dimensionless parameter determined by the mass of the particle, the correlation 
length and the average intensity of the field. Depending on the value of this parameter, the system 
exhibits different regimes, characterized by the localization properties of the eigenfunctions. We 
calculate the corresponding density of states using the statistical properties of the speckle potential. 
We find good agreement with the results of numerical simulations. 
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I. INTRODUCTION 

Speckles are high-contrast fine-scale granular patterns 
occurring whenever radiation is scattered from a surface 
characterized by some roughness on the scale of the ra- 
diation wavelength. Originally discovered for laser light 
in the early 1960s, the speckle phenomenon plays an im- 
portant role not only in optics, but also in other fields, 
where it is used for several applications, including, for 
example, radar and ultrasound medical imagery [l|. In 
recent years, optical speckles have been employed in com- 
bination with cold atoms Q and especially Bose-Einstein 
condensates (BECs) @ in order to investigate the behav- 
ior of matter waves in the presence of disordered poten- 
tials Many interesting features of BECs in one- 
dimensional (ID) speckle potentials have been addressed 
from both the experimental and the theoretical sides. 
These include classical localization and fragmentation ef- 
fects, frequency shifts and damping of collective excita- 
tions, and inhibition of transport properties [H-[l6j], and 
have culminated with the observation of Anderson local- 
ization for a BEC fl2l - fl4l ]. Also, the superfiuid- insulator 
transition fl7l - l24j and the transport of coherent matter 
waves have been recently investigated from the theoret- 
ical point of view for speckles in higher dimensions [25l — 

& 

Despite this intense research activity, the properties of 
the single particle spectrum of the speckle potential have 
been addressed only partially, even in the ID case. In 
[To[ | it has been argued that for low energies the den- 
sity of states (DOS) of blue-detuned repulsive speckles is 
characterized by a Lifshitz tail [29| , whereas the presence 
- for higher energies - of an effective mobility edge, has 
been discussed in [l2| . In addition, though there is a vast 
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literature about random ID systems [2^, [301, most of the 
general theorems apply strictly to the case of uncorre- 
cted disorder, whereas the speckles, being a correlated 
disordered potential, deserve an explicit treatment. 

In this article we discuss the properties of DOS for 
a ID speckle potential, by comparing analytical predic- 
tions, based on the statistical properties of the speckles, 
with explicit calculations for the spectrum of numeri- 
cally generated speckle patterns. We consider both blue- 
detuned (repulsive) and red-detuned (attractive) cases. 
The speckles are characterized by their intensity in and 
a correlation length £, which represents the length scale 
over which the modification of the potential occurs. 
Therefore, for dimensional reasons, the single-particle 
properties are determined by a single dimensionless pa- 
rameter s = ±2m^ 2 /o/fi 2 , with m being the mass of the 
particle. In our conventions, the plus and the minus signs 
refer to the blue-detuned and red-detuned speckle, re- 
spectively. 

In the presence of disorder, the asymptotic behavior 
of the DOS near the edges of the spectrum is dominated 
by rare large fluctuations of the potential. For the red- 
detuned potential, the low-energy tail of the DOS origi- 
nates from the deep wells at negative energies. The latter 
is closely related to the distribution of intensity maxima 
in the speckle patterns. In the case of a blue-detuned 
speckle the spectrum is bounded from below. The DOS 
near the edge results from the large regions with very 
small intensity. For both cases we analyze the behavior 
of the DOS throughout the s range, from the semiclassi- 
cal limit |s| — > oo, down to the quantum regime \s\ < 1, 
where smoothing effects become important. 

The article is organized as follows. In Sees. [TT| and IIIII 
we introduce the speckle potential and its statistical 
properties. Then in Sec. II VI we discuss the DOS for the 
attractive red-detuned speckles, for different regimes of s. 
The case of the blue-detuned repulsive speckles is consid- 
ered in Sec. [Vj Our results arc summarized in Sec. IVII 
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II. THE SPECKLE POTENTIAL 

The electric field £{x) created by a laser speckle on 
a ID Euclidean line at the image plane is, to a very 
good approximation, a realization of a complex Gaussian 
variable. In general, £(x) can be taken as statistically 
isotropic in the complex £ plane and statistically homo- 
geneous and isotropic on the x line. Thus, (£(x)) and 
(£(x)£(y)) vanish while the autocorrelation function is 
given by Q, HH 



(£*(x)£(y)) = G(x -y)=I sine [D(x - y)] . 



(1) 



where sinc(t) = sin7rf/(7rf) and D is the aperture width. 
The probability distribution of intensity I(x) = \£ (x)\ 2 
across a speckle pattern follows a negative-exponential 
(or Rayleigh) distribution 0, HH , 



P[7] = exp[-//J ]// 



(2) 



where Iq = (I) is the mean intensity while the most prob- 
able intensity is zero. The intensity autocorrelation func- 
tion is given by (I(x)I(0)) = 7 2 [l + sinc 2 (Dx)] . The 
(auto) correlation length £ of the speckle potential is de- 
fined through the equation sinc 2 (D£/2) = 1/2, which 
gives the width at the half maximum of the autocorre- 
lation function. It is related to the aperture width by 
£ = 0.88/D. 

When the speckle electric field is sinned on a sample 
of atoms, this results in a disordered potential felt by the 
atoms, that is proportional to the local intensity of the 
speckles, V(x) = al(x) with I(x) = \£(x)\ 2 . In general, 
if the wavelength of the light is far detuned from the 
atomic resonance, the constant a is proportional to the 
inverse of the detuning A, and can be either positive or 
negative [32| . This corresponds to a disordered potential 
bounded from below and composed by a scries of barriers 
(A > 0, blue detuning) or bounded from above and made 
of potential wells (A < 0, red detuning) . In what follows, 
we include |a| in redefinition of Iq, which is positive as 
well as I(x). Then we use a = +1 and a = —1 for the 
blue-detuned and red-detuned speckles, respectively. 

The quantum properties of a particle of mass m in the 
speckle potential are determined by the solutions of the 
Schrodinger equation, 



ip(x) = Ei(j{x), 



(3) 



with ip( x ) being the particle wave function. Note that 
there are two relevant length scales in the problem [33j . 
The first one is the correlation length £ of the disorder 
which introduces the corresponding energy scale, 



E/: = h 2 /2m£ 2 



(4) 



Additionally to this length scale, the particle of mass m 
moving in the random potential with typical strength Iq 
introduces a new length scale B defined by 

I = h 2 /2mB 2 . 




Figure 1. (Color online) The typical ID (repulsive) blue- 
detuned speckle profile. The intensity is given in units of 
Iq and the coordinate in units of £. 



By comparing Iq with the dimensionless parameter s 
can be written as 



a BY 



(5) 



For large values of |s| the disorder is strongly correlated 
while for \s\ — > the limit of uncorrelated disorder is ap- 
proached. Expressing the intensity of the speckle poten- 
tial in units of Iq so that (/) = 1 and the length in units 
of £, we can rewrite the Hamiltonian in a dimensionless 
form as H = — V| + sl(x) with a single parameter s. 

Once the electric field correlation function ([1} is 
known, the statistical properties of the intensity are com- 
pletely determined. This is discussed in the next section, 
where we consider the statistics of extrema of this corre- 
lated and non-Gaussian potential, relevant for the motion 
of a particle in a random landscape [34|. By using this 
description, we derive the DOS and compare the analyti- 
cal predictions with the spectrum computed numerically 
for a randomly generated speckle pattern. 

To generate the speckle potential profile we use the 
mathematical counterpart of the mechanism that gives 
rise to optical speckles. The constructive or destruc- 
tive interference of randomly phased elementary compo- 
nents of the radiation field caused by the roughness of 
the scattering surface occurs also in the discrete Fourier 
transform of a sample function of a random process 
This effect can be exploited to generate numeri- 
cally a speckle pattern with the statistical properties dis- 
cussed earlier 0, [3l| . Then, the spectrum of the system 
can be computed numerically by mapping the station- 
ary Schrodinger equation Hip = Etp on a discretized grid 
[35| . All the results presented in this article have been 
obtained for a system of length L = 600£ and are aver- 
aged over a number of realizations ranging from 10 2 to 
10 3 . 
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III. STATISTICAL ANALYSIS OF THE 
SPECKLE POTENTIAL 

The typical shape of the potential profile created by 
a blue-detuned speckle pattern is shown in Fig. [TJ The 
potential profile can be seen as series of peaks (maxima) 
statistically distributed around the typical value of the 
order of Iq. The peaks are separated by valleys of typical 
width £ with minima situated most probably near the 
lower-boundary of the potential. For the red-detuned 
case the picture is analogous, with the roles of maxima 
and minima exchanged. In other words, the statistics of 
the potential minima for the red-detuned case is that of 
the maxima of the blue-detuned one. 

We are interested in the density of minima n min (/') 
and the density of maxima n max (I'), at specified in- 
tensity I', in the spatial distribution I(x) of intensities 
along the image line. Using the method of Weinrib and 
Halperin [36j , the number of minima or maxima can be 
calculated from the integrated joint probability, 



N, 



mm, max 



dx 5 (I x ) I xx 9 (±I XX 



where the plus (minus) sign refers to minima (maxima) 
and we have used the conventions 81/ dx = I x and 
d 2 I/dx 2 = I xx . It follows that the density of minima 
(maxima) per unit length and intensity can be written 
as 



0.20 
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o . 05 
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(/') = (S(I X )S(I'-I) \I XX \Q(±I XX )). (6) 



Figure 2. (Color online) Density p(7/7o) of minima (top) and 
maxima (bottom) for a blue-detuned speckle (in log scale in 
the insets). The solid lines represent the analytical prediction 
of Eq. I|12p . while the dots have been obtained by averaging 
over 100 numerical speckle realizations. 



In order to evaluate the statistical average of Eq. ([6]) , we 
need to calculate the joint probability of I, I x and I xx (all 
at the same point) starting from the joint probability of 
the Gaussian random variables £, £ x , and £ xx . Because 
£ (x) is statistically homogeneous along the x coordinate, 
the density in Eq. ([5]) will be independent of position. 

The statistics of a set of Gaussian variables is com- 
pletely determined by the two-point correlations G(x~y) 
between them. The nonzero correlators evaluated at 
the same point arc calculated by considering the deriva- 
tives of G(x - y). One finds (\£ | 2 ) = G(0), (£ xx £*) = 
-(\S X \ 2 ) = G xx (0), (\£ xx \ 2 ) = G xxxx (0) with 



G(0) = io, 

Gxx(u) = -i /~e, 

G xxxx (0) = Io/(k £ 4 ), 



(7) 
(8) 
(9) 



where we have introduced the length £ = \/3/(irD) for 
convenience. The dimensionlcss parameter 



G xx {$) 



G(0)G xxxx (0) 



(10) 



is restricted in one dimension to the range < k < 1. It 
depends only on the the shape of the correlation function 
of the electric field but not on the length scale £. For the 



correlator of the laser speckle given by Eq. ([T]), we have 
K = 5/9. 

On dimensional grounds, it is clear that the density of 
minima (maxima) of Eq. © scales as 



^mnijinaxV 



(I, ^ K) = l/(JoO Pmi„,max (J/I 0) k) . (11) 



Evaluating the statistical average in Eq. ((5]) and rewriting 
in terms of dimensionless intensity measured in units of 
Iq, we obtain 



-'mm, max 



I,k) = 



1 



4^ 3 / 2 /1/2^ K -1 - 1 

/"CO 

de x f(i,e x ) di xx g (i,e x J xx ), (12) 

Jo 



-iel (^_ 3 + K - 1 )/( K - 1 -i) 



where / = exp 

g = I xx exp {-I xx [(l xx /4 fj T (0 X - l)] / (k^ 1 - l)}, 

and the plus (minus) sign refers to the minima (maxima). 
The details of the derivation are given in the Appendix. 
The densities (fT2"j) are plotted in Fig. [5] The agreement 
with the numerical simulation is excellent. 

The asymptotic behavior of the integrals in Eq. (fl~2")) 
can be also evaluated analytically. This is useful because 
the distribution of the deep negative-energy wells of the 
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Figure 3. (Color online) The normalized distributions of dis- 
tance measured in units of the autocorrelation length £ be- 
tween neighboring maxima (dotted line) and between neigh- 
boring minima (solid line). 



red-detuned speckle potential is related to the large I 
asymptotic tail of the maxima distribution of the blue- 
detuned speckle. This latter behaves for k = 5/9 as 



Pmax(-0 



2?r 



1 + 21) erf[l] 



V7 




(13) 



Near the boundary of the blue-detuned speckle potential 
the density of minima diverges as 



,(/) « 1/2VttJ. 



(14) 



However, in the thermodynamic limit, the DOS of a par- 
ticle in disordered potential bounded from below is ex- 
pected to vanish at the boundary. Hence, the correlation 
between the DOS and the number of minima has to be 
rather weak. 

More likely, the DOS near the boundary of the spec- 
trum is associated with the distribution of distance 
between neighboring minima P m i n (L) and neighboring 
maxima P max (L). It is known that the presence of Lif- 
shitz tails in the DOS for a bounded-from-below potential 
is related to the existence of large regions free of disorder 
for the binary distribution or large regions with negli- 
gible potential for continuous distributions. In the case 
of blue-detunded speckle potential, the zero intensity is 
the most probable so that one would expect that such 
regions arc not so rare. Some idea about their statis- 
tical properties can be gained through the distributions 
of distances between minima and maxima, which can be 
computed numerically. The normalized distributions of 
the distance between maxima and between minima are 
shown in Fig. [3] The average distance for the both distri- 
bution is (L) = 1.78 £. This value can be easily obtained 
by integrating the density in Eq. (|T2"j) over all the in- 
tensities and by assuming that the points of minima (or 



maxima) are homogeneously distributed. The right tails 
of both distributions can be well fitted in the window 
L/£ G [3,5.5] by simple exponentials: 



P max (L) « 2.081e- 1 - 933 ( L /«), 
P min (L) « 3.385e- 2 - 036 ( L /«). 



(15) 
(16) 



The exponentially small probability of existence of a re- 
gion of size L with relatively small intensity accounts for 
the Lifshitz tail in the DOS of the blue detuned speckle 
potential. 



IV. DOS: RED-DETUNED SPECKLE 

We consider first the DOS for a particle in the red- 
detuned speckle potential, which corresponds to s < 0. 
We discuss separately three different regimes, ranging 
from the semiclassical limit, for s — > — oo, down to the 
quantum regime, when |s| < 1. 

By varying the parameter s the characteristic exten- 
sion of the localized states differs considerably. In Fig. |4] 
we present the lowest-lying cigenstates for three different 
intensities of the same speckle profile, corresponding to 
s = —0.1, —1, —100. This picture shows that for shallow 
potentials (the quantum regime) the eigenstates extend 
over several potential wells. The lowest-lying eigenstates 
localize inside deep wells. However, the ground state 
is not necessarily localized in the deepest well since the 
width of the well plays also a crucial role. As the speckle 
intensity is increased, the eigenstate width shrinks and 
their positions may change. Eventually, for large enough 
|s|, the eigenstates tend to stack up as bound states in- 
side isolated wells, with the ground state being the lowest 
eigenstate of the deepest well (semiclassical regime). 



A. Semiclassical regime 

In the limit s —> — oo the low-energy tail of the DOS is 
determined by deep isolated states occupying single wells 
of typical size £. Therefore, for E < 0, the DOS can be 
calculated in the semiclassical approximation j37| 



,{E) = T I* ME -I) e-l 7 '/ /0 dl, 

1 J-QO 



(17) 



where vq(z) = yfm / '^flithsfz is the density of state of 
the continuum in one dimension. Integrating over the 
intensity in Eq. (|17[) we find 



2ith 2 h 



(18) 



The DOS in the negative-energy tail predicted by 
Eq. (|18[) is in excellent agreement with the curve obtained 
in the numerical simulation as shown in Fig. [5] Note that 
this result differs from the DOS of a long-range correlated 
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Figure 4. (Color online) Probability density distribution of 
the lowest lying eigenstates of a red- detuned speckle pattern 
(bottom panel) for three different intensities corresponding to 
s = -0.1, -1, -100. Only a window of length 150£ of the 
whole system is shown. The figure shows the ground state 
(blue solid line) and several first excited states (black dashed 
line). Note that, as the speckle intensity changes, the position 
of the eigenstates may change and some may move in or out 
from the selected window. 



Gaussian potential that, sufficientl y d eep in the tail, has 
the form v{E) ~ exp(— const|i?| 2 ) [29| . 

The semiclassical result of Eq. (|18|) can be also un- 
derstood in terms of the density of minima of the red- 
detuned speckle potential. Since the latter is the blue- 
detuned speckle potential taken with negative sign, the 
distribution of minima is nothing but n max , namely, the 
distribution of maxima of Eq. @ studied in Sec. Hill 
We assume that each deep state occupies an isolated sin- 
gle well associated with the corresponding minimum of 
the potential. This approximation is not true in gen- 
eral but is valid for deep enough states. The poten- 
tial inside of the well can be completely described by 
derivatives taken at the bottom of the well, denoted 
by {/} := {I,I X = 0,I xx ,I xxx ,...}. Each well is char- 
acterized by M [{I}} bound states with energies E n , 
n = 1, . . . , M , which can be found from the solution of 
the Schrodinger equation for the particle in the potential 



I(x) = -{I+±I xx z 



It 

6 x 



+...). 



The DOS can be then calculated by summing over all 
levels in each well and over all wells as follows 



M[{i}] 



v{E) 



J ^ — n 



n^ n {{I})5(E-E n [{I}\). (19) 
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Figure 5. (Color online) Large negative-energy DOS ex- 
pressed in units of (£^o) -1 in the presence of a red-detuned 
speckle potential with s = —2000 as a function of the rescaled 
energy E/Iq in linear scale and in logarithmic scale (inset). 
The dots represent the numerical data (averaged over 100 re- 
alizations), while the solid line indicates the analytical semi- 
classical curve (1181). 



For large £ the typical wells which dominate the DOS at 
large negative energies can be considered within a har- 
monic approximation by setting I(x) ~ —I — \l xx x 2 , 
with I xx < so that the resulting spectrum corre- 
sponds to a quantum oscillator with the frequency u> = 
yj— I xx /m. The summation over n can be taken to infin- 
ity because in the limit \s\ — > oo the levels of the harmonic 
oscillator approach a continuum. Therefore, for E < 
we can approximate (fT9)) as 



fHO (E) i 

{I x )5[-I + hw\ 



n=0 



1. 



e (-ix»)).(20) 



We have computed the statistical average and carried 
out the summation over n in Eq. (|20[) by summing over 
large numbers of states. The obtained curve is found 
to be exactly on the top of the semiclassical curve of 
Eq. ([T5]) shown in Fig. [S] Corrections beyond the har- 
monic approximation and/or correlations between differ- 
ent points could be taken into account, including higher- 
order derivatives of I{x) at the bottoms of wells. 



B. Intermediate regime 

Upon decreasing |s| the quasiclassical approximation 
becomes inadequate, as shown in Fig. [6] In general, at 
sufficiently large negative energies E <C Eg, the DOS is 
determined by the deep wells of size A < £ no matter how 
short the correlation length £ is [38]. At higher energies 
E > Eg, the wave function of the localized state is spread 
over more than one typical fluctuation of the disorder 
potential, that is, A > £. This picture is supported by the 
numerical analysis of the spatial extension of the localized 
states in the low-energy tail. As one can see from Fig. @] 
the wave functions for s = —0.1 corresponding to the 
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Figure 6. (Color online) The dots represent the numerical 
data in linear scale and in logarithmic scale (inset) for the 
low-energy DOS tails in units of (Bio) 1 in a red-detuned 
speckle potential as a function of E/Iq for different values 
s < 0: s — —2000 (red squares), s = —1000 (blue circles), 
s = — 1 (open circles) and s = —0.1 (green triangles). The 
black solid line refers to the analytical semiclassical result (|18p 
valid at large \s\. 




Figure 7. (Color online) The dots represent the DOS for a 
red-detuned speckle potential with s = —1 in logarithmic 
scale and in linear scale (inset). The solid line describes 
the modified semiclassical approximation of Eqs. (|22|) - (I24I ) 
with /3 = 0.5 valid in the low-energy tail. At positive ener- 
gies, the DOS approaches asymptotically the ID continuum 
vq(E) = y/m/-Kh\/2 E, indicated by the dashed line. 



energy range shown in Fig. [5] spread over a region A ^> £. 
The same analysis for s = — 1 reveals that the typical 
extension of the wave function for the low-energy states 
is of the same order as £. Note that, since the condition 
E <C corresponds to E/Iq <C — X/| s[ , it is not possible 
to see the transition to the deep localized levels at E <C 
E% from the data of Fig. [5] for |s| < 1. 

Here we study the low-energy tail of the DOS for mod- 
erate values of |s| ~ 1. This value occurs typically in the 
experiments with cold atoms in optical speckles [32| . In 
this regime the semiclassical approximation breaks down 
due to quantum-mechanical effects despite that the in- 
teraction between different wells remains negligible. The 
decrease of \s\ results in a depiction of the DOS seen in 
Fig. [SJ This effect can simply be explained by the reduc- 
tion of the number of levels in each potential well. In the 
harmonic approximation (|19p . one does indeed find that 
decreasing £ at fixed B leads to an increase of the dis- 
tance between levels. Alternatively, reducing Iq at fixed 
£ causes a shrinking of the interlevel distance which is 
counterbalanced by a larger decreasing of the wells depth. 
Thus, the decrease of s always corresponds to a reduction 
of the number of bound states in each well. 

In this case we cannot rely on the infinite-sum ap- 
proximation. Taking into account only a ground state 
in each well severely underestimates the DOS. Keeping 
any constant finite upper limit M[{/}] in the sum over 
levels remains inadequate due to a quite broad distribu- 
tion of the number of bound states in different wells. A 
different approach has to be considered. One method 
is to extend the semiclassical approximation given by 
Eq. (jTT)) in order to include quantum effects. To this 
end it is useful to introduce the ID cumulated DOS [39[ 
V(E, V) = (V2m/irh) ^E~V which is related to the 
total number of states in a well of size A created at the 



minimum V m in = — ^max of the potential V(x) = —I(x), 

N(E, V min ) = (v^to/tiTi) - V~ A. 

The condition to have at least one state requires that 
this number is larger than 1, that is, > E > V m - m + 
H 2 /2m\ 2 . This can be incorporated in the cumulative 
DOS by shifting the energy zero [39j : 



V(E, V min ) 



V m 



r, 2 



2m\ 2 ' 



(21) 



Changing from the cumulated DOS to the DOS we arrive 
at 



v \s\~l(E, A) 



VZirh J- 



P[-V] 



E- V - 



dV, 



(22) 



where P[— I] is the normalized probability distribution 
function of minima of the red-detuned speckle potential, 
so that P[I] is the similar distribution of maxima of the 
blue-detuned speckle: 



P[I] 



"■max 



f 

Jo 



(23) 



Moreover, because we know that the extension of the 
wave function of these states does not exceed the typi- 
cal size of the potential wells, we can take the energy of 
the localized state proportional to the zero-point energy 
inside the well, that is, 



h 2 /2m\ 2 - /3 \E\, 



(24) 



where the parameter /3 cannot be determined exactly 
within this method. However the result is not extremely 
sensitive to the choice of f3 in the limits < j3 < 1. 
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Choosing the value f3 ~ 0.5 the modified semiclassical 
method based on Eqs. (f^^j) and (f2"B")) reproduces very well 
the data for s = -1, as shown in Fig. [7] We also tried 
to fit the computed DOS by a stretched exponential. We 
found that it cannot be done well in the whole window of 
numerical data. However, in smaller windows it can be 
fitted by a stretch exponential with the exponent varying 
between 1 and 1.5 and weakly depending on /3. 



C. Quantum limit 

When the correlation length of the disorder £ is re- 
duced to \s\ < 1, the wave function of the localized states 
in the low-energy tail — l/\s\ < E/Iq < is spread over a 
distance A larger than the typical distance between two 
minima (see Fig. [4}. Therefore, the probability of con- 
structing a local wave function is not directly related to 
the probability distribution of the bare potential P[I] or 
that of its extrema P[I]. The interplay between different 
wells has to be taken into account as well. This can be 
done by considering the probability P\ [I] of the potential 
integrated over the characteristic extension \{E) of the 
state of energy E . 

For finite X(E) > £ we can divide the region occupied 
by the wave function into n = A(P)/£ parts, which are 
approximately uncorrelated. Assuming that the corre- 
lation inside each part is strong enough we can approx- 
imately replace the distribution of potential integrated 
over the correlation length P(. [I] by the bare distribution 
©. Then the distribution integrated over the width of 
wave function P\[I] can be computed using the charac- 
teristic function method. The characteristic function is 



defined as fit) 
gives 



(e M ) so that f x (t) = [f s (t/n)] n . This 



(n) 



-nl/Io 



(25) 



with n = A(P)/£. It can be shown that for very large A 
the integrated exponential distribution (|25]l approaches 
asymptotically a Gaussian law: 



Px[I] 



A»l 



Va 



exp 



Hi - Iq 



(26) 



By reversing the argument, we can say that for a given 
energy there is an optimal well of size A created by the 
average potential V, which may bind a state whose zero- 
point energy should not be larger than h/2m\ 2 . There- 
fore, the negative low-energy tail of the DOS can be cal- 
culated with a variational principle from [391 ] 



s|<l 



(E,X) 



Px\-V\ 



E-V 



2mA- 



dV. 



(27) 



By making the variation with respect to A, for large neg- 
ative energies \E\ > I , we find that the DOS of Eq. (gTJ) 
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Figure 8. (Color online) The dots show the low-energy 
DOS tail computed numerically for a particle in red-detuned 
speckle potential with s = —0.1. The black solid line corre- 
sponds to the variational result of Eq. (|27[1 with Pa [I] given 
by Eq. (|25|) . The inset shows the number of wells occupied 
by the states as a function of energy. The green squares and 
red circles are computed for the first few states in a particu- 
lar speckle sample using the attenuation factor of 10 and 100, 
respectively. The solid line is the prediction of the variational 
calculations. 



goes asymptotically as 



v \a\<.i{ E ) ~ CX P 



8V3H\E\ 3 ^ 
9V2 v/mfig 



(28) 



The limit s — > of very short disorder correlation length 
corresponds thus to the white-noise potential 



(V(x)V(x')) Px = k 2 5(x - x'), 



where Kq 



Iq£ has to remain finite as £ — >• 0. There- 
fore, we see that the asymptotic DOS (|28l) reproduces the 
typical Lifshitz tail law for unbounded delta-correlated 
potential 



(E) - exp 



'_2\/3 f\E[ 
9 Ul 



3/2N 



(29) 



1 /3 

where E^ = (mn^/h 2 ) /2 is nothing but the ID Larkin 
energy [H,E|. 

In Fig. [6] we show some numerical results for the sys- 
tems with values of s down to s = —0.1. For such a 
small value of |s|, the typical wave- functions of the states 
in the low-energy tail are spread over a distance a few 
times larger than the typical distance between neighbor- 
ing wells. These states are not yet sufficiently extended in 
order to reach for the distribution of Pa [Pj the asymptotic 
Gaussian probability shape of Eq. (f2l>| . The probability 
distribution Pa[^] given by Eq. (|25|) lies somewhere be- 
tween the exponential law ^) and the Gaussian curve of 
Eq. (f2"rJ|) . To compute the DOS for moderately small s, 
we substitute Eq. (|25p into Eq. (|27[) and minimize it with 



8 




0.4 r 



0.3 



a 0.2 



0.1 



0.0 



■ • 

■ ■• 

• ■ 



0.30 
0.25 
0.20 
0.15 
0.10 
0.05 



0.00 1 

0.0 



4 



Figure 10. (Color online) The one particle DOS for a blue- 
detuned speckle potential computed numerically for s = 0.1 
(greed circles) and s = 1 (magenta squares). The black solid 
line is the DOS of free particle. Inset. The best simultaneous 
fit of the Lifshitz tail of both DOS by Eq. (|36l) with the same 
c = 0.8141 and A = 0.4098. 



Small s limit 



Figure 9. (Color online) Probability density distribution of 
the lowest lying eigenstates of a blue-detuned speckle pattern 
(in the bottom panel), for three different intensities with s = 
0.1, 1, 100. Only a window of length 150£ of the whole system 
is shown. The figure shows the ground state (blue solid line) 
and several first excited states (black dashed line). 



respect to the number of occupied wells n = A/£. This 
gives the result which interpolates between the semiclas- 
sical result of Eq. ^ and the Lifshitz tail (J2SJ). The 
corresponding DOS computed for s = —0.1 is shown in 
Fig. [5] together with numerical data. The variational 
prediction for the number of occupied states is shown in 
the inset. The dots in the inset depict the number of oc- 
cupied wells for the first few states in a particular speckle 
sample. To estimate this number from the numerical data 
we calculate the width of each wave function using the 
criterion of the amplitude attenuation by factors of 10 
and 100. We find a good agreement with the results of 
the variational calculation. 



V. DOS: BLUE-DETUNED SPECKLE 

The random potential created by a blue-detuned laser 
speckle is bounded from below so that the DOS vanishes 
at energies lower than the natural boundary Eq = 0. 
Analogously to the red-detuned problem, the localized 
region of the spectrum is controlled by the parameter s. 
The lowest-lying eigenstates for three different intensities 
of the same speckle profile, corresponding to s = 0.1, 1, 
and 100 are depicted in Fig. [SJ This picture shows that 
with decreasing s the wave functions extend over larger 
numbers of potential wells. 



The DOS computed numerically for s = 0.1 and s = 1 
are shown in Fig. [TU] At large energies the DOS ap- 
proaches the limit of free particle while at small energies 
it exhibits the Lifshitz's tail behavior. It is well known 
that the DOS near the bottom of a bounded-from-below 
disorder potential is controlled by the existence of very 
large regions "free" of disorder potential. The proba- 
bility of their appearance in the speckle pattern can be 
estimated as follows. 

We assume the existence of a bound state with energy 
E localized inside this region. In general, this region has 
the average intensity that is not exactly zero but well 
below the energy level E. Since we are interested in the 
asymptotic low energy behavior, we can fix the largest 
intensity to E <C Iq. However, very small intensity corre- 
sponds also to a vanishing electric field. The probability 
to have a region of size L with the electric field in the 
interval between £(x) and £ (x) + d£{x) is given by 

dW[£{x),d£{x)] 



exp 



£*{x)G- 1 {x ~ x')£{x')dxdx' 



x(detG)- 1 TT S^M, (30) 
JiOL] ™° 

where ao is the UV cutoff and the inverse correlator is 
defined by 

dx"G(x - x")G-\x" - x') = 8{x - x'). (31) 

The Fourier transform of the electric field correlator ([T]) 
reads 



G(k)=I /DQ(irD-\k\), 



(32) 



9 



where Q(x) is the Heaviside step function. In what fol- 
lows, we use D ~ l/£ absorbing the numerical factor of 
0.88 into the overall coefficient. For small s the eigen- 
states spread over distance L much larger than the typ- 
ical well width £. In this regime det G can be estimated 
using Eq. (|32|) as follows: 



dctG 



-ir/L 

n 



lot 



n 



lot = 



2L/i 



(33) 



k— — 7r/^ k—7v/L 



Substituting this into Eq. (|30[) and putting £(x) — in 
the exponential gives the probability for the region of size 
L to have the electric field in the interval between and 
dS(x) = E l ' 2 e ie , where 9 is arbitrary. Integrating out 9 
and summing over x we obtain 



W[0,E^ 2 } ~ (I 0~ 2L/ H^E/a ) 2L/ao . 



(34) 



Since we are interested in the leading exponential behav- 
ior, we can replace in Eq. (|34[l ap by £ and arrive at 



W[0,E^ 2 ] 



exp 



2L fl 



(35) 



The region being almost free of disorder has a bound state 
with energy E = CqH 2 / '(2mL 2 ) , where c is a constant of 



order 1. Using the latter in Eq. (|35]l to express L in 
terms of E 1 , we obtain the Lifshitz's tail for the energy 
distribution as 



v(E) sa A exp 



-cos 



-1/2 



In 



(36) 



The logarithmic corrections to the Lifshitz tail of the 
same form as in Eq. (|36j) was found for the uncorrelatcd 
disorder by mapping the problem to the random walk of 
a classical particle in a medium with traps [42| and using 
a variational approach j43|. We fit using Eq. ([36]) the 
lower energy tail of the DOS computed numerically for 
s = 0.1 and s = 1. The corresponding curves are shown 
in the inset of Fig. [TO] We find that taking into account 
the logarithmic corrections is necessary in order to fit the 
data in the energy interval that is accessible numerically. 



B. Large s limit 

The DOS computed numerically for s = 100, 1000, and 
2000 are shown in Fig. [TT] The DOS exhibits a quite flat 
crossover region between the narrow Lifshitz's tail and 
the large energy asymptotic behavior corresponding to 
the free particle. For general s the Lifshitz tail of the 
DOS is expected to have the form 



v(E)~exp -coy/h/ET(aM(h/E)) 



where !F(s,t) is a universal function depending only on 
the form of G(x). Fitting the DOS computed numerically 
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Figure 11. (Color online) The DOS of a particle in a blue- 
detuned speckle potential are computed numerically for s = 
2000 (red diamond), s = 1000 (blue circles) and s = 100 
(green squares). The black solid line refers to the analytical 
semiclassical result (|37[1 valid at large s and E/Io ^ l/\/s- 
Inset. The best simultaneous fit of the Lifshitz tail of three 
DOS by v(E) w Aexp -c s~ 1/4 yfh]Ehi(I / ' E)\ with the 
same c = 0.6375 and A = 0.3425. 



(see inset of Fig. [TTJ) we find that for s £ [f 00, 2000] this 
function can be approximated by .F(s, t) fts s -1 / 4 t. Thus, 
the correlations in the speckle potentials growing with s 
renormalizc the effective exponent of s in the exponential 
of the Lifshitz tail. 

To compute the DOS in the crossover regime, we use 
the semiclassical approximation and obtain 

v s ^oo(E) =i / ME- I) e- ///o dl 

2 Jo 



27Tfr 2 / 



-B//o erfi 



(37) 



where erfi(x) is the imaginary error function. The 
DOS (|37"1) goes asymptotic as ~ vq(E) for large positive 
energies. Near the boundary of the spectrum E — > it 
behaves as 



The comparison between Eq. (|37|) and the numerical data 
data in Fig. [TTJ shows that there is a range of energies 
where the two results agree quite well. 



C. Effective mobility edge 

The possible existence of the effective mobility edge 
in ID systems with certain correlated disorder has been 
considered in several works [HI, HI]. It is expected, for ex- 
ample, for the potentials with correlations having a finite 
support in Fourier space. In particular as follows from 
Eq. (|32[) the speckle potential has no Fourier components 
with \k\ > tt / D. As a result, in the Born approximation 
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there is no backscattcring for particles with momentum 
larger than Hk > Htt/D that defines the effective mo- 
bility edge. The localization length above this edge is 
determined by corrections to the Born approximations 
and larger by one order of magnitude. In our unities the 
effective mobility edge is expected at 



E 



M 



/ [1 + (0.88tt) 2 /s] 



(38) 



which gives Em/Iq ~ 8.64 for ,s = 1 and Eyi/h ~ 1.008 
for s = 1000. At these energies the DOS is already very 
close to the DOS of a free particle. One cannot see any 
signature of the effective mobility edge like it happens 
with the DOS and the real mobility edge in higher di- 
mensions. To probe the effective mobility edge, one has 
to study the asymptotic behavior of the wave functions 
at this energy |12| which is beyond of the scope of this 
article. 



VI. SUMMARY 

In this work we have studied the DOS of a quantum 
particle in a laser-speckle potential using the statistics of 
the speckle potential profile and numerical simulations. 
The DOS in the speckle potential can be characterized by 
a single parameter s, which depends on the correlation 
length of the disorder, the average value of the laser in- 
tensity and the mass of the particles. We have considered 
both cases of red-detuned and blue-detuned potentials for 
which the random potential is bounded from above and 
from below, respectively. 

For the negative-energy region of DOS in the red- 
detuned speckle there arc three different regimes: semi- 
classical (|s| ^ 1), intermediate (|s| ~ 1), and quan- 
tum (|s| < 1). For strongly correlated disorder, \s\ ^> 1, 
the tail of the DOS is related to the statistical distri- 
bution of the deep minima of the disorder potential. In 
the quantum limit with \s\ < 1, the DOS is controlled 
by smoothing effects leading to the Lifshitz tail of un- 
corrected Gaussian potential. In the crossover regime, 
\s\ ~ 1. the DOS can be described within the scmiclassi- 
cal approach by including quantum corrections. 

In the case of a blue-detuned speckle potential, the 
DOS near the bottom of the spectrum is given by the 
Lifshitz tail related to the existence of large regions with 
negligible intensity. We have found logarithmic correc- 
tions to the celebrated Lifshitz result due to the con- 
tinuous nature of the distribution of the potential. The 
relevance of these corrections for typical setup parame- 
ters in current experiments is supported by our numerical 
results. 
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Appendix 

Because E x has nonzero correlations only with itself, the 
total joint probability in Eq. © can be written as 

p(E,E x ,E xx )=p(E x )p(E,E xx ). (A.l) 
In rescaled unities of Iq and £ we can immediately write 

p(E x ) = {i?(E* x E x ))-V*e-v$a* = Tr-V^-I". 

(A.2) 

The joint probability for the two correlated variables E 
and E xx can be written using matrix representation as 

P (E,E XX ) = (^ 2 dct M)-V AtM ~ lA , (A.3) 

with At = {E*,E* X ) and 

M= f (EE*) (EE* XX ) \ _ ( 1 -1 \ 
iVi " ^ (E XX E*) (E XX E* XX ) J { -1 k- 1 J ' 

which has an inverse matrix that reads 
M _1 = — - — ( 1 

K- 1 ~ 1 V 1 1 

Equation (|A.3[) can be rewritten as 



(A.4) 



(A.5) 



p{E,E xx ) = 4- e -^(l g l a "- 1 +^-l a + a ^[ g ' g »]), 

hi 1 

Using the following variable transformation E(x) = 
\/l(x) and setting I x = 0, we obtain 

E x =i9 x I^ 2 e ie ^ (A.6) 

( Ix 



TT - I xx 

xx \ 2/V 



--J^je^ + i^/y, (A.7) 



Therefore, the joint probability (jA.lj) can be written as 
p(E, E X ,E XX ) = J p{0, 9 X ,9 XX ,I, I x = 0, I xx ), (A.8) 
where 



p(9, 9 X , 9 XX , I, I x — 0, I xx ) — ——: — - 



(A.9) 



+16^+1^-2161 
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and the Jacobian is Jo = 2~ 3 . Integration over 9 yields 
a factor 2tt. We can also analytically integrate out / and 
9 XX , which is a Gaussian integral. As a result, we arrive 
at Eq. Unj). 
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